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Abstract 

There has recently been interest in relating properties of matrices 
drawn at random from the classical compact groups to statistical charac- 
teristics of number-theoretical //-functions. One example is the relation- 
ship conjectured to hold between the value distributions of the characteris- 
tic polynomials of such matrices and value distributions within families of 
L-functions. These connections are here extended to non-classical groups. 
We focus on an explicit example: the exceptional Lie group Gi. The value 
distributions for characteristic polynomials associated with the 7- and 14- 
dimensional representations of G2, defined with respect to the uniform 
invariant (Haar) measure, are calculated using two of the Macdonald con- 
stant term identities. A one parameter family of L-functions over a finite 
field is described whose value distribution in the limit as the size of the 
finite field grows is related to that of the characteristic polynomials asso- 
ciated with the 7-dimensional representation of G2. The random matrix 
calculations extend to all exceptional Lie groups. 
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1 Introduction 



Most work on the connection between random matrix theory and L- 
functions has concentrated on random matrices chosen from ensembles 
related to the classical compact groups. Montgomery [19], Rudnick and 
Sarnak [23], and Bogomolny and Keating [1, 2] calculated the correlation 
functions of the zeros of the Riemann zeta-function, scaled to have unit 
mean spacing, in the limit as T, the extent of the averaging range up the 
critical line, tends to infinity. Their results suggest that these correlation 
functions coincide with ones relating to the eigenvalues of unitary matrices 
in the limit as the matrix size, N, tends to infinity. In the latter case the 
average is defined with respect to the uniform invariant (Haar) measure 
on the unitary group U(N); that is, with respect to the Circular Unitary 
Ensemble (CUE) of Random Matrix Theory (RMT). There is extensive 
numerical evidence in support of this connection [20], which is expected 
to extend to the zeros of any given principal L-function. Katz and Sar- 
nak [15] conjectured that the distributions of low-lying zeros in families 
of L-functions are the same as those of the eigenvalues of matrices from 
the various classical compact groups (e.g. the orthogonal group O(N) and 
the symplectic group USp(2N), as well as U(N)), the particular group in 
question being determined by the symmetry of the family. This is also 
supported by numerical evidence [22]. 

It was suggested by Keating and Snaith [16] that the value distribu- 
tion of a given principle //-function on its critical line coincides, in the 
limit as T — > co, with the value distribution of the characteristic poly- 
nomials of random unitary matrices, defined again by an average with 
respect to Haar measure for U(N), in the limit as N — > oo. The random 
matrix value distribution was calculated in [16] by expressing the group 
average in terms of an integral over the matrix eigenvalues, using a for- 
mula for the measure due to Weyl [28], and then relating the resulting 
(iV-dimensional) integral to one evaluated by Selberg. This idea was later 
extended in line with the Katz-Sarnak philosophy to relate the value dis- 
tribution within a given family of L-functions at the centre of the critical 
strip to the value distribution of the characteristic polynomials associated 
with elements of the appropriate classical compact group in the N — > oo 
limit [17, 3, 4]. Again, the random matrix calculations were performed 
using Weyl's integration formula and the Selberg integral. One interest- 
ing feature of these calculations is that in all cases the logarithm of the 
characteristic polynomial, normalized appropriately (by logiV), satisfies 
a central limit theorem in the limit N — > co. This is in agreement with a 
theorem of Selberg which states that the logarithm of the Riemann zeta 
function, normalized appropriately (by log log T), also satisfies a central 
limit theorem in the limit T — > oo. For further related developments see 
[10, 11, 5]. 

It is in this context that we now ask whether there is a connection 
between L-functions and random matrices from the non-classical groups. 
A particularly interesting class of these groups, closely related to the clas- 
sical groups, is that of the exceptional Lie groups. Our purpose in this 
note is to point out that a number of the key constructions which serve 
to provide the link with random matrix theory in the classical case have 
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analogues for the exceptional Lie groups. We illustrate this by computing 
the moments and value distribution of the characteristic polynomials of 
matrices associated with the 7- and 14-dimensional representations of one 
particular exceptional Lie group, G2. The methods employed are again the 
appropriate Weyl integration formula and generalizations of the Selberg 
integral conjectured by Macdonald (and known as Macdonald constant 
term identities) [18], proven for G2 by Zeilberger [29] and Habsieger [9] 
and by Opdam [21] in the general case. These methods extend to all of the 
other exceptional Lie groups. We then go on to describe a one-parameter 
family of //-functions over a function field, whose value distribution co- 
incides with that of the characteristic polynomials associated with the 
7-dimensional representation of G2 in the limit as the size of the finite 
field grows (this was proved by Katz [14]). The link with finite fields is 
natural, because N is fixed for the exceptional groups and the L- functions 
in question (whose zeros correspond to eigenvalues) are polynomials. 

This note is organized as follows. In Section 2 we review properties 
of G2 necessary for our random matrix calculations. These calculations 
are performed in Section 3. The L-functions associated with G2 arc con- 
structed in Section 4. In Section 5 we conclude with a brief discussion of 
the generalization to the other exceptional Lie groups. 

2 Preliminaries about G2 

As pointed out in the Introduction, the exceptional Lie groups are closely 
related to the classical matrix groups. One particularly natural way of 
seeing this relationship is via their Lie algebras. An important class of 
Lie algebras, because they form building blocks of more general algebras, 
is the class of complex semi-simple Lie algebras. This class allows a com- 
plete categorization and is elegantly summarized in the possible Dynkin 
diagrams which encapsulate the allowed root systems; the root systems 
describe the structure constants of the Lie algebra (standard references 
for this material include [7, 8, 6]). The result of the analysis is that the 
structure of possible root systems is highly constrained. Indeed the only 
possibilities fall into four infinite families, a„, b n , c n and d n plus five 
exceptional cases, g2, ft, ee, ei and e$. Each complex Lie algebra has 
a compact real form and this real form is the Lie algebra of a compact 
group; the compact real forms of and d n are the Lie algebras 

of SU(n + 1), SO(2n + 1), Sp(2n) and SO(2n) respectively. The five 
exceptional cases are the Lie algebras of the groups G2, -F4, Ee, E7 and 
Es- What is significant for our purposes is the many of the constructions 
which exist for the classical groups have analogues for the exceptional 
groups. 

Here, for concreteness, we will mainly focus on the smallest exceptional 
case, namely G2, for which the Lie algebra is the compact real form of g2- 
The group G2 is the automorphism group of the octonions, and has an 
embedding into SO (7). 

The group is 14-dimensional and has rank 2. The six positive roots 
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may be taken to be 




The complete set of roots is R — {±o?i}, i — 1...6. The set of short roots 
is Rs = {±Qi, ±CK3, ±04}; the set of long roots is Rl = {±«2, ±05, ±ae}. 
The Weyl group W of G2 is a dihedral group with 12 elements. 

Let T be a maximal torus of G2, which is isomorphic to a product of 
two circles S 1 x S 1 . Every element of G2 is conjugate in G2 to an element 
of T, which is unique up to conjugation by the Weyl group. 

Weyl's integration formula reads: If d/j, 1Ilv (g) is the Haar probability 
measure on G2, dt is the Haar probability measure on T and F is a 
continuous function on G2, invariant under conjugation, then 

F(g)d^ inv (g) = -L ^ F(t)\A(t)\ 2 dt (2) 

where 

A(t) = (deta)t' (4) = t 4 I] t 1 - ( 3 ) 

(the equality is Weyl's denominator formula), where 8 is half the sum of 
the positive roots: 

8 = i a = 5qi + 3q 2 . (4) 

a>0 

and a (5) means the vector obtained from 8 by the Weyl group element 
a. If we parametrize a particular maximal torus by t = (ti,t2), then t™ is 
an expression of the form t™' 61 ^'" 2 where ei and e2 are (two-component) 
vectors; the vectors ei and €2 and range of the parameters t\ and t2 
depend on the choice of maximal torus. We derive our results below 
without needing to make an explicit choice for this torus. 

For any pair of integers [m , U2\ there is an irreducible representation 
P[m,n 2 ] which has highest weight \ ni ,n 2 ] = niuii + n2LJ2 where wi = 04 
and u>2 = O-q are the fundamental weights. 

The character xa of the representation evaluated at the group element 
t G T is 

X x(t)=Tr[px(t)] = ^d^, (5) 

where fj, are the weights of the representation, d M is the multiplicity of the 
weight fi. 
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The character xa(^) an d dimension d\ of the representation p\ are 
given by Weyl's formulae: 

XAW = E^(deta)t-C) ; (6) 

* - n °>° (A : j)Q . 

The sums in (6) are over elements a in the Weyl group. With these 
conventions, the orthogonality relation for characters is 



J |A(t)j 2 X[ni,n 2 ](*)X[mi,m 2 ] (*) ^ = <5 ni , mi (5„ 2 , m2 . 



(8) 



We will be particularly interested in the fundamental representations 
[1,0] (induced from the embedding of G2 into SO (7)), and [0,1] (the 
adjoint representation) , which have characters 

x ( i,o]W = 1+ £ * a ; (9) 

aen s 

Xl o,i]W = 2+£t"; (10) 

aS-R 

and dimensions 

d[l,0] = 7; d[0,l] = 14. (11) 



3 Characteristic polynomials 

We will focus on the characteristic polynomials 

Z(U P , 6) := det(l - U p e- t0 ) (12) 

of (unitary) matrices U p coming from a given representation p of the group. 
The group elements which these matrices represent can be thought of as 
being chosen randomly from the group with respect to the uniform invari- 
ant (Haar) measure. We will calculate explicit expressions for the averages 
(over the group, with respect to Haar measure) of \Z\ a for complex num- 
bers s (see [16, 17] for analogous calculation relating to U(N), O(N), and 
USp(2N)). 

First let us consider the modulus of Z. We wish to calculate 

< \Z(U P ,8)\ S > G2 = J \Z(U p (g),9)\ s d^ inv (g) 

= J \dct(l-U p (g)e- lB )\ s dp inv (g). (13) 

Since the integrand is a class function, this integral reduces to an integral 
over the maximal torus T: 

< \Z{U p ,e)\ S >G 2 = y 2 j T \A(t)\ 2 \dct(l-U p (t)e-* e )\ s dt. (14) 
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3.1 The seven-dimensional representation 

This representation is induced by the embedding of G2 as a subgroup of 
50(7). From (9), we can calculate that 

Z(f/ [1 , 0] ^) = det(l-[7 [1 , 0] (t)e- <9 ) = (l-e- ie ) [] (1 - t a e~ ie ). (15) 

aGR S 

We see that Z has a zero at 9 = 0, as is the case for the characteristic 
polynomial of any odd-dimensional orthogonal matrix. Let us define Z as 

Z(U [lfi] ,9) = (1 - e- ie )- 1 Z([7 [1 , 0] , 9). (16) 

We now present a formula for < \Z(U[i }0 ],9)\ s > at 9 — 0. 
We have that 

<\Z(U [m ,0)\ s >G 2 = ± / |A(t)| 2 | [J (l-t a )\ s dt 

= T2 [ i A ^i 2 n (i-*t*. (17) 



but 



and so 



where 



iawi 2 = n c 1 -*")' 

qS-R 



<iz(f/ [1 , 0] ,o)r> G2 =-^ / na-o^, 



(19) 



[1 if a £ Kl 

Consider first the case that s is an integer. The value of the integral 
(19) is then the constant term in the expression 

-LjI(i-0*«. (21) 

The value of this constant term is in turn one of Macdonald's celebrated 
constant term conjectures [18], proved for G2 by Zeilberger [29] and Hab- 
sieger [9] (see Opdam [21] for a uniform proof): 

(3k s + 3k L y.(2ks)\(2k L )\(3k L y. 

12(2fc s +3k L y.(k s + 2k L y.(k s + fci)!(fc 5 )!((fcL)!) 2 ' 1 ' 

where ks (resp. fcz,) is the value of k a for the short (resp. long) roots. 

Thus for the representation [1, 0], ks = s + 1 and kh = 1, and so for s 
a positive integer or zero, 

*\7(TT (3s + 6)!(2s + 2)! 

<|Z(tf [li0] ,0)| >^- (2s + 5)!( s + 3)!( S + 2)!( s + l)!- (23) 
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It follows from Carlson's theorem (see [27]) that then 

s\7(TT n\\'^ T(3s + 7)r(2s + 3) 

< \Z{U [ll0] , 0)1 >G 2 = r(2s + 6)r(s + 4)r(s + 3)r(a + 2) (24) 

for Res > —3/2. To see this, note that one may deduce directly from 
(19) that 2~ 6s < |i((7 [li0] ,0)j s >g 2 is bounded when Res > 0, and from 
Stirling's formula that 

T(3s + 7)r(2s + 3) 

r(2s + 6)r(s + 4)r(s + 3)r(s + 2) 1 ' 



is also bounded in the same half-plane. The function 

(26) 



< \7(u ->r. r(3s + 7)r(2s + 3) 

< \Z{U [lfi] , 0)1 >c 2 - r(2s + 6)r(s + 4)r(s + 3)r(s + 2) 

is therefore regular and bounded in Res > 0, and vanishes when s is a 
non-negative integer. Carlson's theorem therefore implies (24). 

Now that (24) has been established, we may immediately write down 
expressions for the probability density functions associated with the value 
distributions of log \Z(U[i i0 ] , 0)|, 

Pl(x) = ± r r(3^ + 7)r(2iy + 3) 

and \Z(U [1M ,0)\, 

P 2(x) - J- f r(3s + 7)r(2s + 3) 

for any c > 0. One can easily deduce asymptotic properties of the proba- 
bility density functions from these integrals; see, for example [16, 17]. 

We note finally that Z is real and positive at 8 = and so its phase 
there is zero. 

3.2 The fourteen-dimensional representation 

In this case the determinant Z has a double zero at 8 = (corresponding 
to the twice repeated weight 0, see (10)). Thus we define 

Z(U [0A] ,e) = (1 - e- lB Y 2 det(l - U [0A] e~ w ). (29) 

Similar calculations to the previous case show that for s a positive integer 
or zero, < \Z(U[ 0t i] , 0)\ s >a 2 is given by the constant term in 

lj](l-D fc «, (30) 

where in this case, 

k s = k L = s + l. (31) 
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Thus for s a non-negative integer the Macdonald identity quoted above 
implies that 

s\7(TT n\\'^ ( 6s + 6)!(2 S + 2)! 

<|Z(tf [0il] ,0)| >o 2 = 12(5s + 5)!((s + 1)!)3 - ( 32 ) 

Once again, Carlson's theorem may be applied (in this case after mul- 
tiplication by 2 _12s ) to show that 

< imo^W >c 2 = r(6s + 7)r(2S + 3) 3 (33) 

V 10,1)- 71 G 2 12r(5s + g) [p (s + 2 )]3 ^ > 

for Res > —7/6. This can then be used to write down expressions for 
the probability density functions associated with the value distributions 
of log \Z(U[ 0i i], 0)| and |Z(Z7[ 0j i] , 0)|, as in the previous section. 

Also, as in the case of the representation [1,0], the phase of Z is zero. 

3.3 Other representations and 9^0 

For other representations of G2 and for 6 =fc 0, the integrand is not of the 
form 

na-o*-, (34) 

and, as far as we are aware, closed form expressions for these integrals are 
not known. 



4 Value distribution problems over func- 
tion fields 

Our purpose now is to outline the analogy between number fields and 
function fields over a finite field in the context of the value distribution 
of zeta- and L-functions for these cases and the predicted behaviour in 
terms of Random Matrix Theory. 

It was proved by Selberg that the logarithm of the Riemann zeta func- 
tion on the critical line has a Gaussian value distribution [24], and the 
same is true for all L-functions [26] under suitable assumptions. Selberg 
also investigated the "g-analogue" of these results for the value distribu- 
tion of the family of Dirichlet //-functions [25] at a point on the critical 
line and there too obtained a Gaussian value distribution. Precisely, for 
q prime we have q — 2 primitive characters \ modulo q, and for fixed t 
consider the q — 2 numbers 

argL{\ +it,x) 



\J \ log log q 



(X varies over all primitive/nonprincipal characters modulo g). Then as 
q — > 00, these numbers are distributed as a standard Gaussian. 

Our purpose in this section is two-fold: first it is to point out that 
there are corresponding results for various families of i-functions over 



8 



function fields, with the role of the Gaussian being replaced by various 
distributions from Random Matrix Theory; and second to construct an 
example for which the appropriate random-matrix distribution is the G2 
result calculated above. 

Likewise, there are similar results for the moments of the L-functions, 
which in the number field setting are mostly conjectural [16, 17, 4], but 
in the function field setting can sometimes be proved. 



4.1 Zeta functions 

Let A; be a finite field of cardinality q and X/k a (smooth, geometrically 
connected, proper) curve defined over k. The zeta function of X/k is given 
by the series 

Z(X;T)=eMj2 N " — )> 
z — ' n 

n — \ 

where N n — #X(k n ) is the number of points of X over the field k„, the 
extension of k of degree n. The series is absolutely convergent for \T\ < 1. 

Trivial example: take X = P , the projective line. The number of 
points of P 1 over the finite field k„ is + 1 = q n + 1 and so 

Z(P 1 ;T) = 



(1-«T)(1-T)- 
This zeta function has an Euler product 

Z(X;T) = JI(l-T dcgp )-\ |T|<1, 
p 

where p runs over all closed points of X. (In the example of P 1 , the closed 
points p correspond to irreducible monic polynomials p(x) G k[x] with the 
addition of the "point at infinity".) 

It turns out that Z(X; T) is a rational function of T, of the form 

Z(X;T)- P ^ 



(l-T)(l-qT) 



with P(X; T) G 1+TZ[T] a monic integer polynomial of degree 2g, g being 
the genus of the curve X, which we can write as P(X;T) = IljLiU ~~ 
otjT). The inverse roots ctj are thus algebraic integers. Further, there is 
a functional equation T 1— » q/T: 

Z(X;-? f ) = q 1 - 3 T 2 - 23 Z(X;T). 

If we set T = q~ s then the functional equation translates into s 1— > 1 — s. 

The "Riemann Hypothesis for curves over a finite field" (proved in the 
general case by A. Weil) is that all the inverse roots Oj of P(X; T) have 
absolute value ^fq, that is as a function of the variable s all zeros are on 
the line Res =1/2. 

What is especially important for our purpose is that the polynomial 
P(X; T) is the characteristic polynomial of a matrix: there is a unique 
conjugacy class Ox £ U Sp(2g) in the unitary symplectic group such that 
P(X; T) = det(J - q 1/2 TQx)- 
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4.2 Families of curves 



Now consider a "family" of curves X/k. In order to study the behaviour 
of P(X; T) as X varies, it suffices to understand the distribution of the 
conjugacy classes Ox- In several cases, it is known that as q — > oo these 
become equidistributed in USp(2g) (with respect to Haar measure). 

For instance, this is the case for the family A4 g of all fc-isomorphism 
classes of (smooth, geometrically connected, proper) curves of given genus 
g [15, Theorem 10.7.15]. 

This allows one to compute arithmetic quantities such as the mo- 
ments of P(X; T) as X varies in A4 g (k) by using the corresponding (non- 
arithmetic) computation in Random Matrix Theory for USp(2g). Thus 
one finds that for T fixed, say q~ 1 T = 1, one has 

.-oo #M g (k) xe ^ (k) j USp(2g) 

The moments of the characteristic polynomial in USp(2g) were com- 
puted in [17] and are given by 

Jus P (2g) fj[ r(i/2 + j)r(i + s + c? + j) 

The probability density functions for the value distributions associated 
with the polynomial and its logarithm may then be written as integrals, 
as in (3.1) and (28) [17]. In the case of the logarithm of the characteristic 
polynomial, the limit distribution when g — > oo is a Gaussian. 



4.3 L-functions attached to exponential sums 

We consider one-variable exponential sums constructed as follows: let k 
be a finite field with q elements as above, f(x) and h(x) G k(x) be rational 
functions, ip a nontrivial additive character of k (e.g. for k — Z/pZ take 
tp(x) = exp(27riaa;/p), 0/a6 Z/pZ), and x a multiplicative character of 
k x . Set 

S(xb, X ; f,h;q) = ^i>(f(x)) X (h(x)), 

X 

the sum running over all x G k which are not poles of /, h and such that 
h(x) 0. For the finite extension k n of degree n of k, we get nontrivial 
characters ip n = rp o Tr kn / k and Xn = tp ° N fcn / fe by composing with the 
trace and norm maps. Correspondingly we get exponential sums for k n 

Sn(X, i>\ f, h) := S(x n ,i>n, f, h\ q"). 

The L-function is defined as 

00 rpn 

L(S, T) = exp(^ S n ( X , 4>; f, h) — ). 

z — ' n 

These have an Euler product decomposition and are rational functions of 
T. 

In many cases of interest to us, it turns out that L(T) is in fact a 
polynomial of the form det(J — q 1,>2 Os) with 0s a unitary matrix. 
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4.4 Gauss sums 



Given a nontrivial additive character xp oik and a nontrivial multiplicative 
character \ °f k x , one defines the Gauss sum <?(XiV') by 

Correspondingly we get Gauss sums for k n 

gn(x,i>) ■= s(x»i«)- 

To compute the L-function 

°° rpn 

L(g( X , 4>),T) = o X p(J2 9n(x, 1>) — ) 

71 

n — 1 

one can use the Hasse-Davenport relations: 

-g n (x,4>) = (-g(x^)T- 

These give 

L(g( X ,i;),T) = l+Tg(x,i>). 
As is well known, \g(x, 4>)\ = ^fq- Thus we may write <?(x, i>) = ^e t9x 
(we omit the dependence on ip which is of a trivial nature). Setting 
s = 1/2 + id I log q,T= q- s = q- 1/2 e- 10 , we find 

L(9(x,V)>T) = i + e i(9 *- fl) . 

The g — 2 angles {# x : x 7^ Xo} are uniformly distributed in [0, 2ir). 
This is easy to see from Deligne's estimate on hyper-Kloosterman sums, 
see [12, section 1.3.3]. Thus the moments of L(g(x,tp),T) and of its 
logarithm, averaged over x and taken as q — > 00, are the same as those 
for the function 1 + e %6 . 

4.5 Kloosterman sums (i) 

These are the sums 

Kl(o,p) = exp — (xi + x 2 ). 

A P 

xiX2~a mod p 

More generally for a finite field k with q elements, take a nontrivial addi- 
tive character rp and a =fc 0, and set 

Kl(a, q) = ^2 tp(xi+X2)- 

X1X2—CL 

This sum is real (replace x \— > —x), and as Weil proved satisfies 

|Kl(o,«)|<2V9. 

The associated L-function is a polynomial of degree 2 

L(Kl(o, q), T) = 1 + Kl(a, q)T + qT 2 . 

It is of the form det(7 - g 1/2 9 a ) with 6 a € 5(7(2). It was shown by Katz 
[13] that as q — > 00, the q — 1 conjugacy classes {6 a : a G fc x } become 
equidistributed in 5(7(2) with respect to Haar measure. This implies we 
can compute the value distribution of L and logL via RMT on SU(2). 
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4.6 Kloosterman sums (ii) 

We next look at an example of exponential sums in several variables: 
hyper-Kloosterman sums are n-variable sums (n > 2) generalizing the 
previous example, given by 

Kl„(a,g) = ^2 ip(xi + x 2 + ■ ■ ■ + x n ). 

x\x^. . .x n — a 

Replacing ip by ip o Tr fcm / fc gives the sum Kl n (a, q m ). 
The associated L-function is defined as 

L a (T) := exp (-1)" £ Kl n (a,q m )L- . 

\ m=l ) 

This L-function is a polynomial of degree n. It was proved by Katz [13] 
that it can be written as det(J — q ( - n ~ 1 ^ 2 TO n (a, g)), with O n (a, q) G K„, 
where K n is the compact group USp(2n), n even, SU(n), n, q odd, SO(n), 
q even, n / 7 odd, and G2, g even , n = 7. Moreover, as a varies through 
all nonzero elements of k, the q— 1 conjugacy classes O n (a, q) oi K n become 
equidistributed there as q — > 00 while keeping the type of K n fixed. For 
instance, taking n = 7 and g = 2 r , r — » 00 gives 2 r — 1 conjugacy classes 
{Or(a, 2 r ) :0/a£ F2>-} which become equidistributed in G2 as r — > 00. 



4.7 An exponential sum associated to G 2 

Let p be a prime, p > 17, k = Z/pZ, X(2) the unique quadratic character 
(Legendre symbol) of k x , and ^ a nontrivial additive character of k, that 
is t/>(a;) = e 2max/p for some a £ k x . Consider for t € fc x the exponential 
sum 

K.T(t)= Y, X(2)(x)i>(x 7 + tx). 

zGfcX 

These sums were studied by Nick Katz and the results below are due to 
him [14]. 

Note that K,T(t) = X(2)(-l)/CT(t) andso/CT(t) is real if X(2)(-l) = 1, 
that is if p — 1 mod 4, and imaginary if X(2){— 1) = —1, i-e. if p = 3 mod 4. 

In view of the transformation properties under complex conjugation, 
we divide the exponential sum K,T{t) by the quadratic Gauss sum g{X(2)) 
to get a real number. Furthermore, there is a (unique) choice of sign 
e p — ±1 so that 1 

g(x(2)) 

is minus the trace of a matrix 6 t £ SO(7): KLTit) — — TrOj. Moreover, 
this matrix turns out to lie in G2. 

The associated L-function is a polynomial of degree 7, which is a char- 
acteristic polynomial of the element Q t of G2 , 

L(KT'(t),T) = det(/-6 t T). 

As t varies in (Z/p/Z) x , these p — 1 conjugacy classes O t become equidis- 
tributed in G2 as p — > 00. Thus the value distribution of L(JCT' (t),T) at 
fixed T is computed by RMT for G 2 . 



x At the time of writing we do not know how to determine t p . 
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5 Other Lie groups 



The random matrix calculations reported here for G2 generalize straight- 
forwardly to the other exceptional Lie groups. In each case one has a 
Weyl integration formula, which allows the moments of the characteristic 
polynomials associated with representations of the group to be written as 
integrals over the Cartan subgroup, and a Macdonald identity, which en- 
ables the integrals to be evaluated as ratios of F-functions. This prompts 
the question as to whether families of finite-field //-functions can be con- 
structed whose value distributions are given by each of the other excep- 
tional Lie groups. 

As a final remark, we note that in [18] Macdonald gives constant term 
formulae for affine root systems (which are related to Kac-Moody al- 
gebras). This suggests the intriguing possibility of extending the ideas 
described in this paper to families of random matrices arising from the 
representations of the associated infinite dimensional groups. 
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